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Abstract. The fractal shape and multi-component nature of the interstellar medium together 
with its vast range of dynamical scales provides one of the great challenges in theoretical and 
numerical astrophysics. Here we will review recent progress in the direct modelling of interstellar 
hydromagnetic turbulence, focusing on the role of energy injection by supernova explosions. The 
implications for dynamo theory will be discussed in the context of the mean-field approach. 

Results obtained with the test field-method are confronted with analytical predictions and 
estimates from quasilinear theory. The simulation results enforce the classical understanding of a 
turbulent Galactic dynamo and, more importantly, yield new quantitative insights. The derived 
scaling relations enable confident global mean-field modelling. 
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1. Interstellar turbulence 

Apart from stars, the baryonic matter within the Galaxy is in the form of an extremely 
dilute, turbulent plasma known as the interstellar medium (ISM). The multitude of 
physical processes within the ISM entails a rich heterogeneous structure (Spitzer IQTSh. 

Approximating radiative processes by a simplified cooling prescription, and restricting 
the computational domain to a local patch, the turb ulent ISM is now routi nely modelled 



by rneans of three-dirnensional fluid simulations fe.g.-lKorpi et al . 1999: M elioli fc de Gouveia Dal Pino 
2004ISlyz et al.lbOOsHjoung fc Mac Lowll2006l: fPib et aLlbOoi T One main focus of these 
simulations has been to obtain filling factors of the different ISM phases and compare 
them to the classical predictions as well as observations (e . g.. |Dettmar 1992). Further 
topics of interest in clude turbulent mixin g (|Balsara fc Kiml(2005j ^. thermodynamic dis- 
tribution functions fMac Low et al.l|2005|), and line-of-sight integrated column densities 
(|de Avillez fc Breitschwerdt 2005b). 



1.1. The small-scale dynamo 
While various simulations (jKorpi et al. ll999Hde Avillez fc Breitschwerdtl20"05aHMac Low et al 



/uuoi l discuss tne mnuence oi magnetic neias on tne i&M mo r pnoio gy, little is saia aoout 
the actual mechanism of field amplification. Balsara et al. (|2004l ) have addressed this 
question by means of unstratified simulations of SNe turbulence. The authors relate 
the growth of small -scale magnetic fields to vorticity production in super nova shocks 
(jBalsara et al.ll200lh . and chaotic field line-stretching (jBalsara fc Kimll2005h . 

The fact that vorticity production by colliding shells is a lmost inevitable in a clumpy 
and highly structured ISM has first been pointed out bv iKorpi et al. ( 19991 ). The is- 
sue has then been i nvestigated for the simplified case of driven expansion waves by 
Mee fc Brandenburg! ([2006^ and, more recently, by Del Sordo fc Br andenburg (th i s vol- 
ume). Considering turbulence driven by non- helical transverse waves, iHaugen et all (|2004l) 







Supernova- driven dynamo 



1 



have shown that the smaU-scale dynamo becomes harder to excite in the super-sonic 
regime, albeit the critical Reynolds number for the onset of dynamo action only seems 
to depend weakly on the Mach number. 

Because the eddy turnover time is short at small scales, a dyn amo based on cha otic 
field line stretching will be fast. This is in-line with observations fsee lBeck et al.|[T996l and 
this volume), which exhibit dominant turbulent fields. Open issues remain with respect 
to the mechanism governing the saturation of the small-scale dynamo. Therefore, it is 
currently unclear whether equipartition field strengths can be obtained by a non-helical 
dynamo alone. Alternatively, the turbulent field might be explained as a "shredded" 
coherent field, i.e., as the by-product of a helical mean- field dynamo. 



2. The large-scale galactic dynamo 

Notwithstanding the above, the presence of the observed coherent f ields on scales 



large r than the outer scale of the interstellar turbulence ('^ 50 pc, see iFletcher et al 



clearly requires the presence of a coherent dynamo. The favoured ca ndidates for 
such a dynamo are, in no particular order: (i) t he kinetic driving by SNe (jKorpi et al 



19991: iGressel et ar 2008b; Gis singer et"alll2009l). fii ) the buoyant cosmic ray-supported 
Parker insta bility (iParker,,199 Hanasz et al.l 2004[l. (m) the magneto-rotational insta- 
bility (MRI. ISellwood fc Balbusl ll999l: iDziourkevitch et all 120041: iNishikori et al.l liool 
Piontek fc OstrikeJ 2007f ), and (iv) gravitational interactions ( Kotarba et al. 2010l ) 3 

The last of these effects is certainly dominant at the early stages of galaxy formation 
and will provide a seed field for subsequent processes. It remains to be shown, however, 
how important external interactions are in the presence of a re alistic feedback from scales 
currently unresolved in cosmological simulations. Moreover, as Hanasz et al. ()2009 |) have 
shown, the cosmic ray (CR) dynamo critically relies on the anisotropy of the CR diffusion 
coefficient. Combined simulations, including both SNe and CRs, will prove whether this 
anisotropy remains effective for strongly tangled turbulent fields. 

The MRI will be impor tant (at least) in the outer regions of galaxies, where the star 
formation activity is low (|Korpi et al.ll2010l ) - even under moderate turbulence, it may 
operate efficiently. The stability criterion for a global isothermal disk of thickness 2H 
leads to the relation 

8VP^t<Cn = (2.1) 



for instability ([Kitchatinov fc Riidigeij|2004|) . Depending on the local rotation frequency, 
i7, turbulent diffusivity r/t, and turbule nt magnetic Prandtl num ber Pmt - which is 
generally expected to be of order unity ( Fromang fc Stonj 2009() - this value may be 
lower than the critical dynamo number for the supernova-driven dynamo. It should be 
worthwhile to address this question within state-of-the-art simulations of the ISM, with 
an Alvfen velocity of the external vertical field of the order lkms~^. 

2.1. The supernova- driven dynamo 

The energy input through SNe into the ISM is tremendous. The cor responding am plitud e 
of the expected dynamo effect has first been estimated by ISokoloff_ fc Sliukuroyi (Il990h. 
Assu ming hydrostatic equilibrium and applying quasi-linea r theory ( Riidiger fc Kitchatinov 
1993 ). these estimates have subsequently been refined bv iFrohlich fc Schultd ( "9961 ) . A 
shortcoming of the approach was the neglect of a possible galactic wind. In general, there 
was a controversy as to whether the turbulence created by the SNe was too vigorous to 



f Also see the respective reviews by Hanasz, Otmianowska-Mazur, and Lesch (this volume). 
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Figure 1. Snapshots at time t = 72Myr of tiie top iialf of the now larger simulation box 
{upper panels), and the disc midplane {lower panels) . The variables shown are: (a) number den- 
sity [cm~^], (b) column density [cm~^], (c) temperature [K], (d) velocity dispersion [kms~^], 
and (e) magnetic field strength [fiG]. The logarithmic grey scales cover ranges [—4.76, 1.01], 
[17.56,21.83], [2.13,7.03], [-0.61,2.64], and [-5.98,-1.18], respectively. 

warrant dynamo action - either because of a too strong wind (carrying the field away), 
or contrary to this, because of a too strong downward pumping (enhancing turbulent 
dissipation near the mid plane) . 

In a series of papers, Ferrierd analytically derived the electromotive force stemming 
from isola ted remnants: the line of work was als o later supported by simulations of single 
remnants ( Kaisig et al.lll993 : Ziegler et al.lll996h . To obtain the net a effect, a convolution 
with an assumed vertical SN distribution was applied. However, the approach suffered 
from a too weak dynamo and highly dominant (upward) pum ping. Only when consid- 
ering the stratified nature of the galactic disc ( Ferrierel 1998 ). the issue was somewhat 
alleviated, although still predicting a strong upward pumping. 

Pioneering semi-global simulations based on "first p rincip les" where performed by 
Korpi. Brandenburg. Shukurov. Tuominen fc Nordlundl ( 1999t ). and it was only for the 
fact of a too low dynamo number that no field amplification was observed in their sim- 
ulations. The first direct simulations exhibiting dynamo action were reported almost a 
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Figure 2. Dynamo coefficients as a function of tlie lieiglit z. The Variables are plotted in dark 
{aRH, . . .) or liglit (q0^, . . .) colours, respectively. Sliaded areas indicate la-fluctuations. In the 
region for \z\ < 2kpc, the results agree well with a previous run applying a much smaller box. 



decade later ( Gressel et al. 2008bl) . The general morphology of such simulations is illus- 
trated in Figure [TJ note the apparent correlation between the density (a) and magnetic 
field amplitude (e) near the midplane, indicating field amplification via compression. In 
contrast, away from the midplane (i.e., in the diffuse medium) the field shows more folded 
structures. 

The fast growth of the dynamo can be understood rigorously in terms of mean-field 
theory, i.e. via a Reynolds-averaged induction equation 



^fB = Vx£: + V X [ uxB - ryVxB ] 
and assuming a standard parametrisation 

£i = OLijBj — fjijSjkidkBi 



hi e {-R,0},fc = 



(2.2) 



(2.3) 



with tensorial coefficients mj and fjij. The determination of these closure parameters 
(Gressel ct al. ,2008a,) hugely b enefited from the development of the so-called test field 
method (jSchrinner et al.ll2005l ). For kinematically forced turbulence, this method allows 
to determine unambiguously all eight tensor components - and, in fact, none of them 
can be ignored (Grcsscl 2009) . Notably, we find a significant positive Radler effect in the 
off-diagonal elements of the diffusivity tensor (see panel 'd' in Figure [5]) . 

The measured a effect and turbulent diffusion are of the expected sign and magnitude 
(see Figure[2l panels 'a' and ' c') . The key finding is that of a m oderate downward pumping 
(panel 'b') as predicted by iRiidiger fc Kitchatinovl (|l993l ). This inward pumping has 
profound implications in compensating the effect of the equally strong upward mean flow 
(not shown), thus yielding ideal conditions for the dynamo. Moreover , a win d may aid 
the shedding of small-scale magnetic helicity as discussed in lSur et al.l (|2007l) . 
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3. Recent advances 

One major concern with the local box approximation was the issue of its limited dimen- 
sion in the horizontal direction. In the low-pressure ambient medium of the galactic halo, 
supernova remnants can easily expand to several hundred parsec in diameter. Because of 
the periodic boundaries and the fini te domain si z e, thi s leads to (spuriously) correlated 
self-interactions. The simulations of iKorpi et al. I (|l999h used Lx ~ Ly ^ 500 pc, which 
was justified for the limited vertical extent of their model. For our standard runs we 
applied = Ly = 800 pc together with a vertical box size of ±2 kpc. 



3.1. Large box simulations 

Even with near kpc horizontal box size, the issue of self-interaction was still seen in 
the far halo region. To eliminate the possibility of an artificial dependence of the results 
on the horizontal box size, we carried out a fiducial simulation run at = 1.6 kpc 
(cf. Figures [T] & [2]) • In the region of overlap, the extended profiles agre e well with th eir 
counterparts derived from the smaller boxes (cf. Fig. 4.5, model 'F4' in lGressell2009l ). 

Beyond ±2 kpc, the velocity dispersion increases significantly, enhancing the turbulent 
diffusion and resulting in noisy profiles for the other coefficients. Note, however, the pro- 
nounced systematic a effect in the region of strong vertical gradients near the midplane 
(panels 'a' and 'b' in Figured]). This clearly supports the paradigm of helicity produ ction 
being the direct result of inhomogeneous turbulence ( Riidiger fc Kitchatinovlll993h . 



3.2. Scaling relations 

Beyond the purpose of mere diagnostics, the profiles shown in Figure [2] serve as a foun- 
dation for global large eddy-simulations solving (|2.2[) . We believe that the semi-global 
approach of vertically stratified local boxes captures the essential physics behind the 
turbulent a effect in the galactic disc. At the same time, we are aware that the local 
geometry implies certain restrictions on the permitted dynamo modes. Identifying global 
symmetries is however necessary for a direct comparison with observations. 

To warrant global mean-field modelling of the a fl dynamo, the scaling of the mea- 
sured closure coefficients with the relevant input parameters has to be obtained. This is 
important because the angular velocity f2, the local shear rate q, the su pernova rate g , 
and the midp l ane de nsity pc are functions of the galactocentric distance ( Ferrierel 2001 ) 



Gressel et all (|2009[ ) already made a first step in this direction by inferring the depen- 



dence on the super nova rate. Corresponding mean-field models have been reported in 



Elstner et al.l ( 20091 ). Further scaling relations with respect to 51 and pc will be reported 



elsewhere. 

First estimates (subject to small number statistics) for a range of 0.1-1 times the 
galactic star formation rate lead to the relations 

Ca^nH, and Co cr-O-5 . (3.1) 

This has the consequence that the dynamo number D = CaCfi ~ il^-^H^a'^ '^ scales 
inversely with the star formation rate. Therefore we would not expect a stronger am- 
plification of the large-scale field for stronger star formation activity. Lastly, the pitch 
angle p ~ a/ Coi/Cq H~'^-^a'^-^^ , in the kinematic regime, only seems to depend 

weakly on the studied parameters. The strongest dependence here is, however, expected 
from the local shear rate q = dlnil/dlni?, suggesting a dedicated parameter survey. 
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4. Conclusions and prospects 

Because current numerical simulations are limited to very moderate Reynolds num- 
bers, it is important to understand how efficient the observed mechanisms remain under 
realistic conditions. To achieve this, it has prov en fruitful to s tudy sirnplified scenarios 
and run multiple parameter sets (see iBrandenb urg fc Subramanian 2 0051 for a compre- 
hensive review). Fortunately, the emerging physical effects are dominated by the outer 
scale of the turbulence, i.e., as soon as a rudimentary scale separation is achieved, the 
turbulent quantities should become independent of the actual micro scale. 

The growing complexity of models challenges the distribution of computing time: in- 
creasing physical realism leaves little margin for the variation of key parameters, let alone 
convergence checks or running multiple representations of a single parameter set. Ded- 
icated studies_j;emAm_mai^^ artificial trends from genuinely physical 
ones ( Hanasz et al. l l2009HGressel et al.ll2009l) . 



The dilemma becomes even more apparent when looking at the recent trend to per- 
forming "resolved" global simulations. For these, convergence checks are a rare exception. 
While an "enhanced" diffusivity may be sufficiently approximated by the numerical trun- 
cation error on the grid scale, the diamagnetic pumping term certainly is not. Yet vertical 
transport has profound impl ications on the emerging dynamo modes and growth rates 



(see, e.g.. lBardou et al.ll200ll ) 



In conclusion, we advocate a strategy that has been applied with great success in the 
design of aircraft, namely the concept of large eddy (or mean-field) simulations. Global 
fluid simulations currently cannot guarantee scale separation for all relevant physical 
scales. To obtain quantitatively correct results, we therefore believe that a sub-grid scale 
model is inevitable. Present local box simulations are a valuable means to provide a 
rigorous framework for the calibration of such a model. 
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